L 2 STABILITY ESTIMATES FOR SHOCK SOLUTIONS OF 
SCALAR CONSERVATION LAWS USING THE RELATIVE 
; ENTROPY METHOD 

O 

^ ' NICHOLAS LLOLIi 

> 

o 

^ ■ Abstract. We consider scalar nonviscous conservation laws with strictly 

' convex flux in one spatial dimension, and we investigate the behavior 

of bounded L 2 perturbations of shock wave solutions to the Riemann 
problem using the relative entropy method. We show that up to a time- 
Qh ' dependent translation of the shock, the L 2 norm of a perturbed solution 

relative to the shock wave is bounded above by the L 2 norm of the initial 
' perturbation. 

Keywords: conservation laws; relative entropy method; riemann problem; 
L 2 stability; shocks. 

AMS Subject Classification: 35L45, 35L65, 35L67. 

o : 

"sf , 1. Introduction 

°^ 

C*~) ' For scalar nonviscous conservation laws with general flux, it is well-known 

from the theory of Kruzkov [20j that the solution operator for the initial 
value problem 

f d t U + d x A(U)=0; 

{ U(x,0) = U°(x), (L1) 

^ ' forms an Z^-contraction semigroup. As a result, the measure in L 1 (R) of 

the difference of any pair of entropy solutions is non-increasing over time. 
In particular, if <f>(x — at) is a shock wave and U° — <p G then 

\\U(;t)-tf>{.-at)\\ L i m < \\U° - <t>\\ LHm (1.2) 



for all t > 0. While Kruzkov's estimate is only valid for scalar equations 
(Temple [25J), global stability estimates for shocks with respect to the L 1 
metric have also been obtained for hyperbolic systems of conservation laws, 
at least for sufficiently small perturbations of suitably weak shock waves. 
One such result is given by Bressan et al. in [5], where the authors establish 
L 1 stability estimates corresponding to an "almost" contractive semigroup 
structure. 

On the other hand, for systems of conservation laws with a convex ex- 
tension, it is well-known that the relative entropy method developed by 
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Dafermos |10| and DiPerna [H] provides L 2 stability estimates for solutions 
away from shocks. As illustrated in [12] . one gets estimates of the form 



\\U(;t) - U{;t)\\LH[-RM) ^ aebt W U ° ~ U \\L2 {[ -R-st,R + st]), (1-3) 



where U and U are weak and strong solutions, respectively, and where a, 
b, and s are constants depending on the initial data Uq and Uq. However, 
simple examples show that this kind of result cannot hold when U is a shock 
wave or more generally when U is only a weak solution of the conservation 
law. Nonetheless, as we will prove, one can expect similar estimates to hold 
up to a suitable translation of the shock. 

Consider the initial value problem (II. ip for a scalar conservation law in one 
space dimension. Our goal is to prove the the following global L 2 stability 
estimate for shocks. 



Theorem 1.1. Let U° £ L°°(R) and assume U° - <j> £ L 2 (R) where 



with Cl > Cr. Further, assume U is the unique entropy solution of 
for a smooth flux function A : R — » R verifying A" > 0. Then there exists a 
Lipschitz continuous function x : [0, oo) — > R and a constant A(||[/°||l°°; (p; A) 
> such that 



for all t > 0, where a is given by the relation o~(Cl — Cr) = A(Cl) — A{Cr). 

In the same spirit as [5], our result can be characterized as an almost 
contractive variation of Kruzkov's estimate (jl.2[) . In fact, we will prove a 
slightly more general result (Theorem 13. 1|) in Section [3] which takes into 
account all strictly convex entropies associated with However, while 

Theorem 13.11 is interesting in its own right, it is important to keep in mind 
that the estimates gained from the relative entropy method are purely of 
type L 2 , regardless of the specific convex entropy used. 

In a related result, Chen et al. [7j have used the relative entropy method to 
obtain stability estimates for shocks in the context of gas dynamics. Specif- 
ically, the authors establish the time-asymptotic stability of Riemann solu- 
tions with arbitrarily large oscillation for the 3x3 system of Euler equations 
in one space dimension. The present work is another attempt at developing a 
stability theory for shocks using relative entropy techniques, beginning with 
a treatment of scalar equations, as proposed by Vasseur in [13]. For further 
applications of the relative entropy method in the context of fluid dynamics 
and kinetic theory, we refer the reader to the papers pQ EJ El [H [191 EH ET] . 




(1.4) 



\\U(-,t) -(!)(■- at -x(t))\\ L 2 (M) < \\U° - ^|| L2(R) 



(1.5) 



and 




(1.6) 
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Let us clarify a few things regarding estimates (|1.5f) and (|1.6p . First, we 
are interested in controlling the relative entropy globally in time. Therefore, 
unlike the formulation of many large-time stability results (cf. |16[ \T7\ \18\ 
[21] I22j). we are forced to take a time- dependent translation of the shock. 
Also, we should point out that fjl .6f) is only intended to show that the shift 
x(t) has a bound proportional to the size of the initial perturbation. The 
estimate captures neither the Lipschitz continuity of x at t = 0, nor the 
expected large-time behavior of x. Indeed, based on our techniques, it is 
reasonable to expect the time-asymptotic convergence of x(t) to the shift, 
xq, considered in the large-time stability analysis of Il'in and Oleinik |17j . 
at least when the initial perturbation is in L 1 (R). 

The proof of Theorem 11.11 relies heavily on the work of Dafermos; in 
particular, on the theory of generalized characteristics and contingent equa- 
tions considered in [9]. Roughly, the idea is to find curves xi,(t) and xji(t), 
initially positioned at the origin, which preserve the quantities 

£ L (t) = f (U(x,t) - C L f dx, and 

/■oo 

£n(t) = / (U(x,t)-C R ) 2 dx. 

Jx R (t) 

Quite surprisingly, this condition leads to the constraint x^(t) < XL^t) for 
all t > 0, so that (|1.5p holds for all functions x(t) satisfying xn(t) — at < 
x(t) < XL{t) — at. The bound on x then follows easily from the local L 1 
stability of entropy solutions. 

The paper is organized as follows. In Section we introduce the entropy 
and relative entropy inequalities associated with the conservation law (jl.lj) . 
and establish some preliminary estimates related to those quantities. In 
Section EJ we describe our method in detail and present the proof of the 
main theorem. Finally, we include in the appendix an existence result for 
differential inclusions arising in the context of conservation laws. 

2. Relative Entropy Estimates 

Consider the scalar conservation law 

8 t U + 8 X A(U) = 0, (2.1) 

with smooth flux A : R — ► R. We say that U G L°°(R x (0, oo)) is an entropy 
solution of (|2.ip if U satisfies 

dtv(U) + d x G(U) < 0, 

in the sense of measures, for all entropy /entropy- flux pairs (r/, G) S [C°°(R)] 2 
verifying rj" > and G' = rj'A'. If U solves f|2. 1 1) in the classical sense, and 
we consider, associated to each convex r/, the relative entropy function 

rj(U | U) = n(U) - r](U) - n'(U) • (U -U), 
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then an entropy solution U verifies additionally the inequality 

d t [ V (U | U)] + d x [G(U) - G{U) - n'(U) ■ (A(U) - A(U))} (2.2) 

< -d t [ V '(U)] -(U-U)- d x [n'(T7)] ■ (A(U) - A(U)), 

for the same entropy/entropy-flux pairs (77, G). 

The idea of the relative entropy method is to use (|2.2p to estimate the 
quantity f n(U \ U) dx in time. When n(U) = U 2 , this corresponds directly 
to estimates in the I? metric. In this paper, we will only consider the case 
U(x,t) = C, for constant states C G R, so that (12.2f) reduces to 

d t r)(U I C) + d x F(U, C) < 0, (2.3) 

where 

F(U, C) = G(U) - G(C) - r,'(C) ■ (A(U) - A(C)) 
u 

(rj'(w) -rj'(C))A'(w)dw, (2.4) 

represents the flux of relative entropy. 

The goal of this section is to establish some preliminary estimates related 
to relative entropy inequality (|2.3p . Before presenting those results, let us 
remind the reader of the following facts. 

Remark. Throughout the paper, we assume that U is an entropy solution of 
(l2Jj) with initial data U° G L°°(R). We assume also that the flux function, 
A : R — > R, is strictly convex and smooth, so that A" o U > a > 0. With 
these assumptions, Oleinik's estimate implies that U(-,t) £ BVi oc (M.) for 
all t > 0. Therefore, the one-sided limits U(x—,t) and U(x+,t) exist and 
verify U(x—,t) > U(x+,t) for all x G R and for all t > 0. Furthermore, the 
trace theorem of Vasseur [26J shows that U(-,t) is countinuous with values 
in Lj oc (R) up to t = 0. 

We begin with the following estimate, which is an adaptation of a key 
lemma in [9] [Lemma 3.2] as it applies to (|2.3p . 

Lemma 2.1. Assume r\ : R — > R is smooth and convex, and letxi '■ [to, T] — > 
R and xr ■ [io>^] ~~ > R be Lipschitz continuous functions such that for all 
t G [to, T], Xft(t) — XL(t) > 5 for a fixed 5 > 0. Then for C G R and for all 
a and b with to < a < b < T, 

C) dx (2.5) 
x L {t)ri(U(x L (t)+,t) I C)dt 

x R {t)r](U{x R (t)-,t) I C)dt, 
for any entropy solution U G L°°(R x (0, 00)) of \2.1\ ). 




x L (b) 



n(U(x,b) I C)dx 



n(U (x, a) 



lx L (a) 

< / F{U{x L {t)+,t),C) 
F(U(x R (t)-,t),C) 
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Proof. For e < |, let ip e (-,t) be an "inner" approximation to the character- 
istic function on [xi(t), x R (t)] (instead of the "outer" approximation taken 
in [9]). Specifically, let 



0, 

\{x-x L (t)) 
1, 

o, 



if x < XL(t); 
if XL(t) < x < x L (t) + e; 
if XL{t) + e < x < x R {t) — e; 
:(*)), if XR(t) - s < x < x R (t); 
if x R (t) < x. 



Also, as in [9] let 



Xe(t) 



0.. 



if t < a; 

if a < t < a + e; 
if a + e < f < b; 
if 6 < t < b + e; 
if 6 + e < i. 



Applying the non-negative test function (p £ (x,t) = tp £ (x,t)xe(t) to (|2.3p and 
taking e —>■ yields (|2.5p . ■ 



Now consider the Riemann data <f> given by (jl.4p . Continuing in the spirit 
of |9j we have the following estimate. 

Lemma 2.2. Let (ft be defined by an d suppose U G L°°(R) satisfies 

~~ ^IIl 2 (K) < +oo. AZso, assume n : R — ► R is smooth and convex, and 
let xl '■ [to,T] — ► R and x# : [toi^l - > R ^ e Lipschitz continuous functions. 
Then for all a and b with to < a < b < T, 

■x L (b) poo "] 

r](U(x, b) | C L ) dx+ j T}(U(x, b) | C R ) dx \ 

-oo Jxn(b) I 

x L (a) roo \ 

7](U(x, a) I C L ) dx+ r](U(x, a) | C R ) dx \ (2.6) 

-oo Jxr(o) J 

F(U(x L (t)-,t),C L ) -x L (t)r](U(x L (t)-,t) | C L )di 
+ / F(U(x R (t)+,t),C R )-x R (t) V (U(x R (t)+,t)\C R )dt, 

J a 

where U is the unique entropy solution of \2.1\) with initial data U°. 
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Proof. Choose R sufficiently large so that —R < xi,(t) — £ and consider 
instead the test function 



0, if x < -2R; 

±(x + 2R), \i-2R<x<-R- 

1, if -R < x < x L (t) - e; 
-\{x-x L (t)), if x L (t) - e < x < x L (t); 



0. 



if xi{t) < x. 



Testing (|2.3p with (p £)R (x, t) = ip £jR (x, t)x £ {t) and taking e — > and R — > oo 
we get 



x £ (6) 



v(U{x,b) | C L )dx 
r6 



xi, (a) 



rj(U(x, a) | Cl) dx. 



< 



lim — 

_R— >oo rt 



-il 



-2il 



F(U(x,t),C L )dx } dt 



F(U(x L (t)-,t),C L )-x L (t)r](U(x L (t)-,t) | C L )d*. 



Since lit/ 



|l 2 (r) < +oo and 



|F(E/(M),C L )| < 



sup |A'(ui)| 
kl<l|t/°l|i-+|Ci| 



(U(x,t)-C L f 



on account of (12, 4j) , the first term on the right-hand side above vanishes. 
The relative entropy on the right is controlled by a similar argument. ■ 



The following corollary is immediate. 

Corollary 2.3. Assume that the hypotheses of Lemma \2.S\ are satisfied, and 
assume additionally that 

(i) x'L,[t) < f{U{xL{t)—,t),C£) for almost every t G [to>^1, and 

(ii) x R (t) > f(U(x R (t)+,t),C R ) for almost every t £ [t ,T], 

where f(U, C) = ^jjpy ■ Then, for all a and b with to < a < b < T, 

x L (b) poo ~\ 

r](U(x, b) | C L ) dx + / r](U(x, b) | C R ) dx \ 

-oo J xji(b) 1 

x L (a) i-oo \ 

tj(U(x, a) | C L ) dx + v(U(x, a) | C R ) dx \ < 0. 

-oo Jx R (a) \ 
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2.1. The Normalized Relative Entropy Flux. Consider the function 
/ : R x R — > R, which we call the normalized relative entropy flux, defined 
by 

f( u > °) = SftttI = I ^ u > C ) A '( W ) dw > ( 2 - 7 ) 

where 



V(U | C) Jc 



The estimates obtained in Section are based on the fact that for strictly 
convex functions A and rj, the function defined by (|2.7h is Lipschitz and 
increasing in the variables U and C, respectively. In order to prove this, 
let us first compute the first order partial derivatives of /. Assume U ^ C. 
Then, 

J£(E7, C) = <p(U, U, C)A'(U) + f ^-(w, U, C)A'(w) dw 

= <p(U,U,C)[A'(U)-f(U,C)} (2.8) 



and 



%{U,C)= f U ^( w ,U,C)A'(w)dw 

J C 



-rf'(C)-{U-C) 
V(U | C) 



(2.9) 



Using these formulas and taking suitable Taylor approximations, one can 
show that / and its gradient have a continuous extension to the line U = C. 
In particular, since 

/ <f(w, U, C) dw = 1, (2.10) 
Jc 

for all U ± C, we define f(C, C) = A'{C) for all Cel. 

As suggested above, we would like to show that on bounded subsets of 
R 2 , (|2.8p and (|2.9p are bounded by positive constants from above and be- 
low, respectively. The latter estimate is delicate and relies on the following 
lemma. 

Lemma 2.4. Let g : [0, 1] — > R and h : [0, 1] — > R 6e continuously dif- 
ferentiable functions such that g'(s) > e 9 > and > £h > /or a// 

s S (0,1). Further suppose Jq g(s)ds = 1. Then, 

i /T _ _ 



g(s)h(s) ds — / h(s) ds > 



o 
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Proof. Since g is increasing and f g(s) ds = 1, there exists a unique sq G 
(0, 1) such that g(so) = 1. Therefore, 



g(s)h(s) ds 



h(s) ds 



(g(s) - g(s ))h(s) ds 



(g{s) - g(so))(h(s) - Hs )) ds 



h'(r) dr 



g'ij) dr 
> £ g £h J (s - s ) 2 ds > e -^- 



ds 



where we used the fact that the last term is minimized at so = \. 
We can now prove our original claim. 



Lemma 2.5. Assume A : 
functions, and let f :lxl 
!! C K x 1, we have 



i and r\ : R — > R are smooth, strictly convex 
be defined by {2. 7|). Then, for bounded sets 



< Lq, and 



H 0<e n <-^ 
where en and are constants depending on Q. 

Proof. First, choose R sufficiently large so that \(U, C)\ < R for all (U, C) £ 
$7, and let ea, £r/, La, and be positive constants such that 

< e A < A"(0 < L A (2.11) 

and 

< e v < rf{0 < L v (2.12) 

for all |£| < R. Taking into account (|2.10p and the fact that A' and rf are 
increasing, the quantity computed in (|2.8p . for U / C, can be controlled in 
the following way. 

df 



< 



nu (U,C)<<p(U,U,C) [A'(U)-A'(C)] 



< 



v(u\c) uc 

Lr,LA(U - Cf _ L V L A 



e v (U - Cf e v 

Since the bound is independent of \U— C\, the estimate holds also for U = C. 
This proves part (i) of the lemma. 
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Next, assume U ^ C and observe that, after a change of variables, equa- 
tion (|2.9h can be rewritten in the form 



-r]"(C) -\U-C\ 



h(s) ds 



where 



and 



v(u | C) 

h(s) = sgn(U - C)A'(C + s(U - C)) 
g(s) = {U- CMC + s(U - C), U, C). 



g(s)h(s) ds 



(2.13) 



Let us check that g and h satisfy the hypotheses of Lemma 12.41 First, 
according to (|2.10p . we have L g(s)ds = 1. Moreover, we compute 



(U - C) 2 r]"(C + s(U - C)) e v (U - C) 



v(u I C) 



2e r . 



> 0. 



and 



h'(s) = \U - C\A"(C + s(U-C))>\U- C\e A > 
Therefore, Lemma 12.41 applies and we deduce from (|2. 13|) that 



d£ (u c) > -v"{c).\u-c\ 
dc {U ' C) - V (U\C) 



2e r , 



\U-C\e A 



12 



12 V L 



£A£ V 



> 0. 



Again, since this bound does not depend on \U — C|, the estimate extends 
to the line U = C. Therefore, we have (ii) and the proof of the lemma is 
complete. ■ 



Note that the proof of Lemma [2.5l in the case rj(U) = tt- is much simpler. 
Indeed, a change of variables in (|2.7p yields 



f(U,C) 



2 f sA'{C + {U -C)s)ds. (2.14) 
Jo 

Therefore, differentiating (|2.14p with respect to U and C and using (|2.1ip 
we get 



and 



2 df 



1 df 



2 r 

< -L A 
3 



< -L A - 
3 
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3. A Relative Entropy Technique for Shocks 

It is well-known that when A is convex and when the initial data (ft, 
given by (11, 4ft , is non-increasing, the unique entropy solution of the Riemann 
problem is the traveling shock wave <ft(x — at), where the shock speed a is 
given by the Rankine-Hugoniot relation. The goal of this section is to show 
that traveling shock waves are stable in the sense of Theorem 11,11 In fact, 
we will show a slightly more general result. 

Theorem 3.1. Let U° G L°°(R) and assume U° - (ft £ L 2 (R) where (ft is 
given by with Cl > Cr. Further, assume U is the unique entropy 

solution of U for a smooth flux function A : R —* R verifying A" > 0. 
Then, for any smooth rj : R — > R verifying rj" > 0, there exists a Lipschitz 
continuous function x : [0, oo) — > R and a constant \(\\U°\\l°° ; (ft', A; rj) > 
such that 

o roo 

n(U(x,t)\(ft(x-at-x(t)))dx< / r)(U°(x) | (ft{x))dx < oo, (3.1) 

XD J — OO 

and 

\x{t)\<\\\U° -(ft\\ L 2 m Vi (3.2) 

for all t > 0, where a is given by the relation o~(Cl — Cr) = A(Cl) — A{Cr). 

In order to show (|3.ip for a strictly convex entropy r], the idea is to 
construct curves xl '■ [0,oo) — > R and xr : [0, oo) — > R, initialized with the 
data xl(0) = Xi?(0) = 0, for which the the total relative entropy 

rx L (t) poo 

8{t)= n(U(x,t)\C L )dx+ n(U{x,t)\C R )dx (3.3) 



J -oo Jx R (t) 

is bounded above by 5(0) for all t > 0. (Note that when U° — (ft G L 2 ( 
and C/° G L°°(R), then 5(0) is finite.) Due to the compressive nature of the 
solution U(x,t), the construction will produce automatically the constraint 
xr{€) < XL{t) for all t > 0. Therefore, (|3.1j) will follow for any function x(t) 
satisfying xn(t) < x(t) + at < xi(t). In particular, this includes the curves 
x(t) = XL^t) — at and x(t) = xn(t) — at. 

To control the total entropy, the idea is to exploit (|2.6p by choosing x'l (t) 
and x'nit) so that the right hand side vanishes. This makes sense at points 
of continuity of U; however, it turns out we do not have as much freedom at 
points where U is discontinuous. We borrow the following lemma from [9]. 

Lemma 3.2. Let x : [to, T] — » R, < to < T < oo be a Lipschitz continuous 
function. Then for almost all t G [toj^"]; 

A(U(x(t)+, t)) - A(U(x(t)-,t)) - x(t) [U(x(t)+, t) - U(x(t)-,t)] = 0. 

The lemma simply asserts that if x(t) moves along a discontinuity of U 
then its derivative must coincide with the shock speed given by the Rankine- 
Hugoniot condition. 
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i(t) 



Next, motivated by the idea of generalized characteristics, we consider a 
curve x(t) solving in the sense of Filippov (see appendix), the differential 
inclusion 

x(t) E [f(U(x(t)+, t),C), f(U(x(t)-,t),C)], (3.4) 

where f(U,C) = ^jf^r ■ I n view of Lemma [3721 we find that (|3.4p is actually 
quite restrictive. 

Proposition 3.3. Let x : [to,T] — ► M be a Filippov solution of [3.$ °n the 

interval [to , T] . Then for almost all t E [to , T] , 

' f(U(x(t)±,t),C), ifU(x(t)-,t) = U(x(t)+,t); 

Ai VmlZ^ ] > ifU( X (t)-,t)>U( X (t) + ,t). 

With these facts in mind, we consider functions xl : [0, oo) — > R and 
xr : [0, oo) — > R, with initial values ccl(0) = x R (0) = , solving 

f x L (t)e[f(U(x L (t)+,t),C L ),f(U(x L (t)-,t),C L )}, 

\ x R (t) e[f(U{x R (t)+,t),C R ),f(U(x R (t)-,t),C R )], '"" J) 

in the sense of Filippov. Existence (and uniqueness beyond t = 0) is guar- 
anteed by Proposition 15.11 (see appendix) and Lemma [2.51 Given (|3.5p . it 
follows immediately from Corollary 12.31 that the total entropy £ (t) in (|3.3p 
is bounded above by £(0). Moreover, since 

r,(U\C)= / v"(0dtdw 



ic Jc 

is non-negative, we easily deduce the following lemma. 

Lemma 3.4. Assume xl ■ [0, oo) — > R and x R : [0, oo) — ► R verify )13. 5|) in 
the sense of Filippov. Further, assume x R (t) < xi(t) for all t > 0. Then for 
any function x : [0, oo) — > R satisfying x R {t) < x(t) < xl{€) for all t > 0, 
we have 

«x(t) roo 

r](U(x, t) | C L ) dx+ I v(U(x, t) \ C R ) dx < £{t) < £(0), 

Jx(t) 

for all t > 0, where £ is defined by A3. 3)) . 

Given Lemma 13.41 it remains to show that when (13. 5p holds and when 
Cl > C R , then x R (t) < x^it) for all t > 0. Since xl and x R coincide at 
t = and they move continuously, the idea is to show that x R cannot pass 
xl', that is, we would like to argue that when x R (t) = XL^t), then in some 
sense x' R (t) < XL(t). This is not precise, of course, since the derivatives may 
be only measurable; however, note that when x R (t) = x^(t) = x and (x,t) 
is a point of continuity of U, then, formally, the monotonicity of / (Lemma 
12.51) implies 

x L (t) = f(U(x,t),C L ) > f(U(x,t),C R )=x R (t). 
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This suggests that xr may only pass ii at a discontinuity of U ; however, as 
we will show, Lemma 12.51 does not allow it. We will argue by contradiction, 
but first let us prove the following lemmas. 

Lemma 3.5. Assume xl '■ \to,T\ — * R and xr : [to,T] — * R verify \3. 5\) in 
the sense of Filippov on the interval [to,T], with Cl > Cr. Further, assume 
that x R (t) - x L (t) >5>0forallt£ [t ,T]. Then for C £ (Cr, C l ) there 
exists A > independent of 5, such that for all a and b with to < a < b < T, 

x R (b) rx R (a) 

7](U(x, b)\C)dx- j r](U(x, a) \C)dx< -\S aJb \\ < 0. (3.6) 

x L (b) Jx L (a) 

where S a ^ = {s £ [a, b] \ x'r(s) - x' L (s) > 0}. 

Proof. Using Lemma \2.1\ for all a and b with to < a < b < T, 

x R (b) rx R (a) 

n{U(x, b) I C)dx- / r)(U{x, a) \ C) dx 



x L (b) Jx L (a) 

< I v(U(x L (t)+, t) | C) [f(U(x L (t)+, t), C) - x L (t)] dt 

V (U(x R (t)-,t) | C) [f{U{x R {t)-,t),C) - x R {t)\ dt (3.7) 
£ / n (U(x L (t)+, t) | C) [f(U(x L (t)+, t), C) - f(U(x L (t)+, t), C L )] dt 

r-b 

n(U(x R (t)-,t)\C) [f(U(x R (t)-,t), C) - f(U(x R (t)-,t), Cr)} dt, 



where we used (13. 5p to get the last inequality. Now, for t G S a fi 

f{U{x L {t)+,t),C L ) < x L (t) < x R (t) < f(U(x R (t)-,t),C R ). 
Therefore, since U £ L°°, we deduce using Lemma [231 that 
f(U(x R (t)-,t),C L ) - f(U(x L (t)+,t),C L ) 

> f{U{x R {t)-,t),C L ) - f(U(x R (t)-,t), C R ) 
c R 9C< 



for t € S a b, where e > is a lower bound on S?. We deduce further, using 



(U(x R (t)-,t),z) dz 

Cn) > 0, (3.8) 

af 

dC- 

and Lemma 12.51 that 

U(x R (t)-,t) - U(x L (t)+,t) > j{C L - Cr) > 0, 

for t G S a ^, where L > is an upper bound on Thus, for t G S a> b, either 
(i) n{U{x R (t)-,t) | C) > K(U(x R (t)-,t) - Cf > « [^(C L - Cr)] \ or 
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(ii) V (U(x L (t)+,t) | C) > K (U(x L (t)+,t) -C) 2 >k [^(C L - C R )\ 



where k > is a lower bound for hrj"(-). Returning to (|3.7p . we get 



x L (b) 



r)(U(x,b) | C)dx 



< 



r,(U(x L (t)+,t) | C) 
-b 



x R (a) 

r](U(x, a) | C) dx 

x L (a) 

C L df 



c 



dC 



(U(x L (t)+,t),z)dz 



dt 



rj(U(x R (t)-,t) | C) 



c 



Cr 



dC 



(U(x R (t)-,t),z) dz 



dt 



< -eK / [ v (U(x L (t)+,t) | C)+ri(U{x R (t)-,t) | G)\ dt 



< -eKk 



2L 



(C L ~ C R ) 



\S n 



(3.9) 



where K = minlC/, — C,C — Cr} (and we used the fact that r](U \ C) is 
non- negative) . This proves (|3.6j) with A = eKk [^(Cl — Cr)] 2 . m 



Lemma 3.6. There exists 1 < k < oo such that for any 5 > 0, if xl : 
[t ,T] —* R and x r : [t ,T] R verify 

' x L (t) £ [f(U(x L (t)+,t),C L ),f(U(x L (t)-,t),C L )}, 

< X R (t) £ [f(U(xR(t) + ,t),C R )J(U(XR(t)-,t),CR)}, 

< x R (t ) - x L (t ) = 6, 

in the sense of Filippov, with Cl > Cr, and if XR{t) — xi{t) > S for all 
t G [to, T], then XR(t) — xi{t) < k5 for all t £ [to,T]. 

Proof. Given t £ (to,T], let St 0t t = {s£ [to,t] \ x'r(s) — x'l(s) > 0}. Then, 

x R (t) - x L (t) = 5+1 [x'r(s) - x L (s)} ds 
J t 

<5+ I [x R (s) - x L (s)] ds. (3.10) 

Since U £ L°°, Lemma [2.51 implies \x'r(s) — x'l(s)\ is bounded by some 
constant Mi. Furthermore, on account of Lemma [3.51 for C £ (Cr,Cl) 
there exists A > such that 

r'x R {t) f-x R {t ) 
0< / n{U(x,t)\C)dx< 7](U(x,t )\C)dx-\S kht \\ 

Jx L (t) Jx L (t ) 

fXR(t ) 

< / M 2 dx-\S kht \X = M 2 5-\S kht \X, 

Jx L (t ) 



14 



N. LEGER 



where again we used U G L°° to bound rj(U \ C) by a constant M2. There- 
fore, |Si 0) t| < <5 (^), and we deduce using A3. lQf) that xr{£) — xi,{t) < k5, 
where k = 1 + ■ 



We can now prove our previous claim. 

Proposition 3.7. Assume xl : [0,oo) — * R and xr : [0, 00) — * R verify 
\3. 5\) in the sense of Filippov with Cl > Cr. Then xu(t) < /or a// 

t > 0. 

Proof. We argue by contradiction. Let d(t) = xji(t) — xi(t) and suppose 
d(T) > for some T > 0. Then for < 5 < d(T) define 

d~ 1 (S) ={0<t<T\ d(t) = 5} 

and let t$ = sup t. Since d is continuous and d(0) = 0, d~ 1 (5) is nonempty 
ted- 1 (5) 

and tg < T. Also, we must have d(t) > 5 for t € [tj,? 1 ], otherwise t$ would 
be larger. Therefore Lemma 13.61 applies with to = ts an d we conclude that 
d(t) < k5 for all t G [ts, T]. In particular, d(T) < k<5 for 5 arbitrarily small. 
This is a contradiction. ■ 



Proof of Theorems 11.11 and 13.11 Since U is the unique entropy solution 
of (jl.ip . the inequality (|2,3p holds in the sense of measures for any strictly 
convex entropy r\. Thus, given xl and xr defined by (|3.5p . the estimates in 
Section [2] are valid, and we deduce from Proposition 13.71 that there exists a 
Lipschitz continuous function x : [0, 00) — > R such that x_r(£) < + at < 
for all t > 0. Also, since U° - <j> € L 2 (R), we have 

L f 00 

r](U°(x) I 0(a;)) dx<^ (U°(x) - 4>{x)f dx < 00, 
—00 2 J— 00 

where > is an upper bound on 77". Therefore, £(0), defined by ()3.3j) . is 
finite and Lemma 13.41 implies 

("OO 

f7([/(z, t) I <f>(x -at- x(t))) dx < £(t) < £(0) = / r](U° (x)\<j>(x)) dx, 



which completes the proof of estimate (13. 1 [> . Also, this gives (ll.5p . for The- 
orem [LTJ in the case T](U) = U 2 . 

Finally, let us show that the function x has a bound proportional to the 
size of XJ — 4> in L 2 (R). Recalling that XR{t) < x(t) + at < xl(£), we have 

\x(t) + at\ < max{\xR(t)\,\x L (t)\} < Lt, 

since xl and xr are Lipschitz (Proposition I5.ip and xl(0) = xr(0) = 0. 
Note that L depends on ||J7°||zo°, (j), rj and A, as the velocities x'l and x'r 
are given by (12. 7p . 
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Next, observe that cj)(x — at — x(t)) — <fi{x — at) has support contained in 
the interval \—{L + \a\)t, (L + \a\)t]. Therefore, 

(C L - C R )\x{t)\ = ||0(- -at- x(t)) ~ 0(- " vt)\\LHB {L+wnt ) 
< ||0(.-(rt-x(t))-^(-,t)|| £l(B(£+k|)t) 
+ ||C/(.,i)-^(.-^)|| I/1(B(i+k|)t) (3.11) 

Then, by the L 1 -stability theory of Kruzkov, the last term is bounded by 
\\U° ~ <t>hi(B (M+L+[irnt ), where M = su V {\A'{w)\; \w\ < \\U Q \\ L ~ + ||0||^}. 
Therefore, proceeding with the estimate (j3.11|) . using Holder's inequality, 
we get 



(C L - C R )\x(t)\ < y/2(L + \a\)t ■ \\<P(- -at- x(t)) ~ U(;t)\\ L2(R) 

+ y/2(M + L + \a\)t- \\U° - <f>\\ip W (3.12) 

Finally, since %(U - <j)) 2 < rj{U \ (f>) < ^f{U - <p) 2 , using (J3l[]) we have 



U{--at-x{t)) -U(; t)\\ L 2 
which together with (|3.12p implies 

\x(t)\ < 



Er. 



Il 2 ( 



1 



Cl — Cr 



1 2L r . 



(L + \a\) + yj2(M + L + \a\ 



)Vi. 



This completes the proof. Note that, by construction, x is actually Lipschitz 
even though this estimate does not show it. 
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5. Appendix: Filippov Solutions and Conservation Laws 

In this section we include an existence result for differential equations 
arising in the context of conservation laws. While this result is an easy 
application of the theory of Filippov [15] . and by no means original, we have 
found no explicit statement of this kind in the literature. 

Definition. A solution of (|5.ip in the sense of Filippov on an interval [to, T) 
is an absolutely continuous function x(t) for which (|5.ip holds for almost 
every t £ [t ,T). 
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Proposition 5.1. Let U be the unique entropy solution of with A" > 
and U° E L°°(R). Also, let g : R —>■ R be continuous and non- decreasing. 
Then for (xo,to) £ R x [0, oo) there exists a Lipschitz continuous function 
x : [in>oo) — ► R, with intial value x(to) = xq, solving 

x(t) G [g(U(x(t)+, t)),g(U(x(t)-,t))}, (5.1) 

in the sense of Filippov (Definition [5|). Furthermore, if to > and g is 
Lipschitz on bounded subsets of R, then the solution x is unique. 

Proof. The existence of solutions follows from [15] [Section 2.7, Theorem 
1] provided the set valued function G(x,t) = \g(U(x+,t)), g(U(x—,t))] is 
upper semicontinuous. Roughly this means that in the limit (x' , t') — > (x, t), 
the sets G(x',t') will be contained in G(x,t). In the present setting, upper 
semi-continuity is an immediate consequence of the continuity of g and the 
structure of solutions to (II. ID detailed in [9] (see also Remark [2]). Also, since 
g o U € L°°, the Filippov solutions are Lipschitz and defined for all £ > to- 
Now assume additionally that g is Lipschitz on bounded subsets of R, 
and let us verify that solutions extend uniquely beyond t = 0. Applying [15] 
[Section 2.10, Theorem 1], it suffices to check that for any T > there exists 
£ € L 1 ({to,T]) such that for any almost every (x,t) and (y,t) in R x [ioi^l 

(x-y)- (g(U(x, t)) - g(U(y, t))) < £(t)\x - y\ 2 . 

The assumptions on g together with Oleinik's well-known decay estimate 
easily imply the statement above. ■ 

Existence results of this type have appeared implicitly in the work of 
Marson and Colombo (see [8] and [23J) on ODEs related to traffic modeling. 
In fact, the uniqueness argument above can be found directly in [8] for the 
case of concave flux functions. The study of contingent equations in the 
context of conservation laws can also be found in the papers [U \9\ lllj. 

References 

[1] F. Berthelin, A. E. Tzavaras, and A. Vasseur. From discrete velocity Boltzmann 
equations to gas dynamics before shocks. To appear in J. Stat. Phys., 2009. 

[2] F. Berthelin and A. Vasseur. From kinetic equations to multidimensional isentropic 
gas dynamics before shocks. SIAM J. Math. Anal, 36(6):1807-1835 (electronic), 2005. 

[3] Y. Brenier. Convergence of the Vlasov-Poisson system to the incompressible Euler 
equations. Comm. Partial Differential Equations, 25(3-4) :737-754, 2000. 

[4] Y. Brenier, R. Natalini, and M. Puel. On a relaxation approximation of the incom- 
pressible Navier-Stokes equations. Proc. Amer. Math. Soc., 132(4):1021-1028 (elec- 
tronic), 2004. 

[5] A. Bressan, T.-P. Liu, and T. Yang. L 1 stability estimates for n x n conservation 

laws. Arch. Ration. Mech. Anal., 149(l):l-22, 1999. 
[6] A. Bressan and W. Shen. Uniqueness for discontinuous ODE and conservation laws. 

Nonlinear Anal, 34(5):637-652, 1998. 
[7] G.-Q. Chen, H. Frid, and Y. Li. Uniqueness and stability of Riemann solutions with 

large oscillation in gas dynamics. Comm. Math. Phys., 228(2) :201-217, 2002. 



L 2 STABILITY ESTIMATES FOR SHOCKS 



17 



[8] R. M. Colombo and A. Marson. A Holder continuous ODE related to traffic flow. 

Proc. Roy. Soc. Edinburgh Sect. A, 133(4):759-772, 2003. 
[9] C. M. Dafermos. Generalized characteristics and the structure of solutions of hyper- 
bolic conservation laws. Indiana Univ. Math. J., 26(6):1097-1119, 1977. 

[10] C. M. Dafermos. The second law of thermodynamics and stability. Arch. Rational 
Mech. Anal, 70(2):167-179, 1979. 

[11] C. M. Dafermos. Generalized characteristics in hyperbolic systems of conservation 
laws. Arch. Rational Mech. Anal, 107(2): 127-155, 1989. 

[12] C. M. Dafermos. Entropy and the stability of classical solutions of hyperbolic systems 
of conservation laws. In Recent mathematical methods in nonlinear wave propagation 
(Montecatini Terme, 1994), volume 1640 of Lecture Notes in Math., pages 48-69. 
Springer, Berlin, 1996. 

[13] C. M. Dafermos and M. Pokorny, editors. Handbook of differential equations: evolu- 
tionary equations. Vol. IV. Elsevier/North-Holland, Amsterdam, 2008. 

[14] R. J. DiPerna. Uniqueness of solutions to hyperbolic conservation laws. Indiana Univ. 
Math. J., 28(1): 137-188, 1979. 

[15] A. F. Filippov. Differential equations with discontinuous righthand sides, volume 18 
of Mathematics and its Applications (Soviet Series). Kluwer Academic Publishers 
Group, Dordrecht, 1988. Translated from the Russian. 

[16] H. Freistiihler and D. Serre. L 1 stability of shock waves in scalar viscous conservation 
laws. Comm. Pure Appl. Math., 51(3):291-301, 1998. 

[17] A. M. II' in and O. A. Oleinik. Behavior of solutions of the Cauchy problem for certain 
quasilinear equations for unbounded increase of the time. Dokl. Akad. Nauk SSSR, 
120:25-28, 1958. 

[18] C. K. R. T. Jones, R. Gardner, and T. Kapitula. Stability of travelling waves for 
nonconvex scalar viscous conservation laws. Comm. Pure Appl. Math., 46(4) :505- 
526, 1993. 

[19] B. Jourdain, C. Le Bris, T. Lelievre, and F. Otto. Long-time asymptotics of a mul- 
tiscale model for polymeric fluid flows. Arch. Ration. Mech. Anal., 181(1):97-148, 
2006. 

[20] S. N. Kruzkov. First order quasilinear equations with several independent variables. 
Mat. Sb. (N.S.), 81 (123): 228-255, 1970. 

[21] T.-P. Liu. Nonlinear stability of shock waves for viscous conservation laws. Mem. 
Amer. Math. Soc, 56(328) :v+108, 1985. 

[22] T.-P. Liu. Pointwise convergence to shock waves for viscous conservation laws. Comm. 
Pure Appl. Math., 50(11):1113-1182, 1997. 

[23] A. Marson. Nonconvex conservation laws and ordinary differential equations. J. Lon- 
don Math. Soc. (2), 69(2):428-440, 2004. 

[24] A. Mellet and A. Vasseur. Asymptotic analysis for a Vlasov-Fokker- 
Planck/compressible Navier-Stokes system of equations. Comm. Math. Phys., 
281(3):573-596, 2008. 

[25] B. Temple. No Li-contractive metrics for systems of conservation laws. Trans. Amer. 

Math. Soc, 288(2) :471-480, 1985. 
[26] A. Vasseur. Strong traces for solutions of multidimensional scalar conservation laws. 

Arch. Ration. Mech. Anal, 160(3): 181-193, 2001. 
[27] H.-T. Yau. Relative entropy and hydrodynamics of Ginzburg-Landau models. Lett. 

Math. Phys., 22(l):63-80, 1991. 

Dept. of Mathematics, 1 University Station C1200, Austin, TX 78712-0257 
E-mail address: nleger@math.utexas.edu 



